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Q Abstract 

^ With their 2010-2011 data set, the LHC experiments have started 

G> their quest to observe the rare decays B^^^ — )■ This study will 

^ provide very sensitive probes of New Physics (NP) effects. NP discovery 

^ potential lies as well in the study of the decay — > K*^iJL^jjL~. Results 

Q and perspectives are presented for studies at the LHC of rare B decays 

^ involving flavor changing neutral currents. 
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1 Introduction 



A powerful way for the LHC experiments to search for new physics is via non Standard 
Model (SM) virtual contributions to rare decays of beauty mesons. The results re- 
ported here for LHCb are based on 37 pb~^ of pp collisions at \fls) = 7 TeV delivered 
by the LHC during 2010. 

2 Search for B^^^ — )► 

Within the SM, the dominant contribution to these modes stems from the Z-penguin 
diagram, while the box diagram is suppressed by a factor of Imw/fntl"^- The Higgs 
annihilation diagram contributes only negligibly (about 1/1000). These are FCNC 
modes which, in addition, are also helicity suppressed. Hence the SM expectation is 
only 3.2 and 0.1 x 10~^ for the Bg and mode respectively [T]. As the error on 
the theoretical expectations is small, these modes are very attractive as a SM test 
bench. They are very sensitive to NP with new scalar or pseudoscalar interactions, 
as well as models with an extended Higgs sector and high tan (3, as their branching 
ratio (BR) goes as (tan/3)^. The current best 95% Confidence Level (C.L.) limits are 
43 and 9.1 x 10~^ for B^ and B^ to respectively [2]; these limits are about a 

factor 10 over the SM expectation. 

LHCb is the only LHC collaboration to date which has presented and published 
a result on the search for these modes, hence its analysis strategy will be described 
in more detail in the following, ii~ pairs are selected on triggered events with 
criteria very close to the ones used for the control samples in order to have efficiency 
corrections as low as possible. The analysis is performed blind. Each candidate 
is given a likelihood to be signal or background-like in a 2D space formed by the 
ensemble of geometrical event variables (GL), and by the invariant mass. To translate 
the number of observed candidates into a BR measurement, normalization samples 
with well measured BRs are used, for which it is necessary to have good knowledge 
of the relative efficiencies. Finally the CLs method [3] in bins of mass and GL is used 
to extract the limits. About 300 events survive the preselection in the mass signal 
region, defined as [Mb^ — 60, Mb^ + 60] MeV/c^. We expect (SM contribution only) 
about 0.3 5° —J- /i and 0.04 B^ — )■ ii'^ [i~ events in 37 pb~^ of data. The power of 
the geometrical discriminant variables is crucial to drastically reduce the rest of the 
background. Given the excellent muon ID, most of the background is actually due 
to real muons from semileptonic decays, while the while peaking background from 
misidentified two-body hadronic B decays is negligible for the integrated luminosity 
at hand. 

Fig. [T] shows the input variables to the GL for signal in red and the main back- 
ground component in black. Such variables, like impact parameters, vertex x^, muon 
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Figure 1: Distribution (LHCb), for signal and background, of the variables used in 
the GL definition. Red: MC fi+fi' signal. Black dashed: MC bb ^ n+fi'X 

background. 
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Figure 2: Distributions (LHCb) of GL input variables for B — )■ hh' (left) and 
K~T[^ (right) decays for data (blue, black squares) and MC simulation (red). Top 
plots show the B (D) impact parameter, middle plots the DOCA of the two tracks, 
bottom plots the isolation variable of one of the tracks. 
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isolation etc, are decorrelated before the construction of the geometrical likelihood. 
Fig. [2] shows some of the variables used in the two-body hadronic B decays B — > hh' 
candidates. In all cases, some differences are observed both in the core and in the tails 
of the distributions. These discrepancies are more clearly visible in the distributions 
of the analogous variables for the high statistics — )■ K~7r^ control sample, also 
shown in the same plot. The optimization and the training of the GL are performed 
on Monte Carlo (MC) simulated data, but the GL output is then calibrated with data, 
thanks to the B — )■ hh' control samples, which have a very similar decay topology, 
hence the small discrepancies between data and simulation do not introduce a bias 
into the analysis. The GL is contructed in a way that it should be fiat for signal and 
peaked at for the background. 
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Figure 3: GL signal (black full squares) and background (blue empty circles) PDFs 
used as input for the calculation of the limit (LHCb). 

Although B-/,hh' decays represent an excellent proxy for the signal mode, they 
are in general triggered in a very different maner to the signal, and this biases their 
kinematical and topological distributions. The solution that has been adopted is to 
to use only events triggered independently of the signal, that is from other particles 
in the event, often those from the other b-hadron. All B — )• hh' channels are fit 
simultaneously, with BR constrained to the PDG values, in bins of GL. The results 
of the calibration, which also takes care of the small differences between data and 
simulation, are shown in Fig. |3]for signal and backfround and it looks very similar to 
what expected from the simulation within uncertainties. 

The GL PDF for the background is obtained from a fit to the mass sidebands in 
GL bins. As expected about 97% of the events are in the first GL bin. 
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Figure 4: Di-muon invariant mass spectrum between 8.75 GeV/c^ and 11.25 GeV/c^ 
at LHCb. The T{1S) resonance has been fitted with a Crystal Ball function while the 
T{2S) and T(35') peaks with Gaussian functions. The background is parametrized 
with an exponential. 



The invariant mass resolution is obtained with two independent methods. The first 
exploits the fact that, according to the Gluckstern parametrization of relative error 
of particle's momenta, the resolution is proportional to the mass of the resonances 
decaying into two muons (after re-weighting the momentum of the muons according 
to those from Bg mesons). Therefore LHCb extracts the invariant mass resolution at 
the mass of the Bg from data by linearly interpolating from the measured resolution of 
charmonium and bottomonium resonances decaying into two muons (i.e. J/ip, ipC^S) 
and T(IS'), T(2S), T{3S)). Fig. |4] shows the T family mass distribution, where all 
resonances are very well separated in LHCb. Another method exploits the B — )• hh' 
modes inclusively. The two results are compatible and are averaged to give 26.7 
MeV/c2. 

For the normalization three different channels have been used: B^^^ — )■ h^h'~, 
B^ — )■ J/ipK~, and Bs — ?■ J/ip(f). Each of these has different efficiencies, different 
dependencies on the Bg and -B^/u hadronization fractions, and different systematic un- 
certainties. For the ratio of B^ and Bg fragmentation fractons, fd/fs, which is needed 
to relate the first two modes to that of the Bg signal decay the HFAG average jl] of 
LEP and Tevatron results have been used ( = 3.71 ± 0.47). 

The numbers of found candidates are in agreement with the expected background 
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level, hence LHCb quotes a limit which is 4.3 x 10^^ for 90% and 5.6 xlO~^ for a 95% 
C.L., slightly lower than the expected average as the distribution of the found events 
are in fact lacking in the very center of the mass distribution, but still compatible 
with the 68% expected band. 

Similar considerations are valid for the — > /i^/i^ mode, where the measured 
limits are 1.2 and 1.5 x 10"^ at the 90 and 95% C.L. 

As for the future, with 200 pb~^, a conservative estimate for this summer, LHCb 
alone could set a 95% limit down to 2 x 10"*^, while with 1 fb~^, expected by the 
end of the year, it will be able to explore the region of BR down to twice the SM 
expectation. A 5 standard deviations discovery is expected for BR values down to 
10~^ with the data foreseen in 2012. 

The analyses at CMS [5] and ATLAS [6] are still blind. They use similar variables 
as LHCb, except for the fact that they use a B isolation variable rather than a muon 
isolation one. The official limit expected for CMS, with the bb cross section currently 
measured, is 2.3 x 10~® for 1 fb~^, but in fact 1.4 x 10~^ according to LHCb modified 
frequentist estimates, hence quite competitive with LHCb. All these expected limits 
have been calculated assuming only the contribution of events from background. 

3 Asymmetries in B'^ — > K*^fi^fi~ decays 

Although the decay i?^ — ?■ K*^fi~^fi~ is much more abundant |3| than that of the 
previous example, it possess high sensitivity to NP through the behaviour of the 
kinematical observables that can be constructed from the final state distributions, as 
the Forward Backward Asymmetry {Aps) in the fifi rest frame, which will be the 
focus of the first analyses with this mode in LHC experiments. What is enticing 
is that at the zero point of this asymmetry the dominant theoretical errors due to 
form factors calculations cancel out. The analyses are performed blind, using control 
samples to limit the dependence on simulation in terms of trigger corrections, selection 
and reconstruction efficiencies, and, especially, acceptance. 

In LHCb the selection is based on Boosted Decision Trees. The clear signal from 
36 pb~^ is shown in Fig. [sj The extrapolated yield for 1 fb~^ of data is about 650 
candidates with a signal over background ratio even larger than for the B-factories. 
All experiments to date see an asymmetry at low opposite, even though statistically 
not significant, to the SM prediction, a possible hint of NP. With 1 fb~^, assuming 
LHCb finds the same central value at low as Belle p], by the end of the year we 
might be faced with an interesting discrepancy of 4 standard deviations with respect 
to the SM, using only S° K*^fi~^fi~ events. 
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Figure 5: LHCb signal for — )■ K*'^fi^fi~ candidates 

4 Conclusions 

LHCb has just published its first 5°/^ — )■ /i^/i^ search results [9J. These results are 
already competitive with the world best limits from Tevatron. LHCb is expected to 
probe the very interesting BR regions very soon. New results from Tevatron, but 
also from CMS and ATLAS, are expected very soon as well. As far as the Aps for 
— )■ K*'^fi~^fi~ decays, LHCb expects to confirm or deny the interesting hints from 
other experiment by the end of this year. 
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